Empirical potential study of phonon transport in graphitic ribbon 
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The thermal properties of graphitic ribbon are investigated based on Brenner's empirical poten- 
tial. The reliability and usefulness of the empirical potential to address the thermal properties of 
covalent-bonded carbon nanostructures are verified through a comparison of phonon dispersion re- 
lations and the density of states of carbon nanotubes with first-principles calculations. The analysis 
reveals unique edge-phonon states that are highly localized at edge carbon atoms of both armchair 
and zigzag ribbons. Applying the phonon dispersion relations to the Landauer formula of phonon 
transport, the quantization and universal features of the low-temperature thermal conductance of 
graphitic ribbon are elucidated, and it is found that the width of the quantization plateau in the 
low-temperature region is inversely proportional to the ribbon width. 

PACS numbers: 44.10.-K, 63.22.+m 



I. INTRODUCTION 

Carbon is unique among all the elements in that it 
can form strong covalent bonds of sp, sp 2 , and sp 3 hy- 
bridizations, depending on the atomic geometry. Owing 
to this strong covalent bond behavior, carbon allotropes 
exhibit remarkable thermodynamic and thermal trans- 
port properties. For example, crystalline diamond has 
the lowest specific heat and highest thermal conductivity 
of all solids. In contrast to bulk crystals, however, it is 
difficult to measure thermal quantities of carbon nanos- 
tructures such as fullerenesi and carbon nanotubes^ as 
it is necessary for measurement to ensure complete ther- 
mal isolation from surrounding materials. The thermal 
transport properties of carbon nanostructures are the key 
to controlling the performance and stability of future 
nanodevices, and it is now widely accepted that ther- 
mal management of nanostructures cannot be avoided as 
the size of devices reaches the nanometer scale. There- 
fore, it has become particularly important to be able to 
understand and predict the thermal transport properties 
of carbon nanostructures through theoretical and com- 
putational analyses. 

The thermal transport properties of carbon nanostruc- 
tures have been widely investigated by molecular dy- 
namics (MD) methods based on the empirical poten- 
tial developed by Brenner^ which is particularly suit- 
able for potential functions of carbon materials. Recent 
MD studies on thermal transport in single-wall carbon 
nanotubes (SWNTs) have predicted that the SWNTs at 
room temperature have relatively high thermal conduc- 
tance, comparable to that of crystalline diamond. 4 ' 5 ' 6,7,8 
High thermal conductance has also been observed experi- 
mentally using micrometer-sized ropes of SWNTs. 9 How- 
ever, the investigation of thermal transport phenomena 
at low temperatures is beyond the capability of classical 
MD methods, as the low-temperature thermal conduc- 
tance is predicted to exhibit quantum effects. Therefore, 
it is desirable to develop a reliable scheme that is ap- 



propriate for the description of low-temperature thermal 
transport in carbon nanostructures. 

The thermal transport properties depend critically on 
the phonon dispersion relations of systems. The phonon 
dispersion relations are obtained by diagonalizing the dy- 
namical matrix containing force constants between any 
two atoms in a system of interest. The force-constant 
matrix is derived from the second derivative of a poten- 
tial function with respect to atom coordinates. This work 
employs the Brenner potential, which has been verified 
to be adequate for the structural study of carbon mate- 
rials including diamond surfaces^ The dynamical matrix 
is derived from the second derivative of the Brenner po- 
tential, leading to the phonon dispersion relations of the 
carbon materials. The phonon dispersion relations for 
SWNTs are calculated first in this paper to verify the 
transferability of the Brenner potential and to demon- 
strate the usefulness of the model. Compared to ab initio 
force-constant methods based on density functional the- 
ory, the empirical model is faster in computation and can 
deal with much larger systems. Most importantly, the 
phonon dispersion relations obtained are in good agree- 
ment with those of ab initio force-constant methods^ii 

The empirical force-constant model is applied in this 
paper to calculation of the phonon transport proper- 
ties of graphitic ribbon as a potential wiring material in 
nanoscale electronic devices as well as carbon nanotube 
applications.* 2 Recently, various forms of graphitic rib- 
bon have been synthesized toward the realization of nan- 
odevicesii&iiiii Polyperinaphthalene is a graphitic rib- 
bon with an armchair-shaped edge (armchair ribbon), 
including 10 carbon atoms in the unit ceMi&H Ow- 
ing to edge effects (finite size), graphitic ribbon ex- 
hibits remarkable physical properties that are not ex- 
pected to appear for carbon nanotubes with closed struc- 
tures. Graphitic ribbon with zigzag-shaped edges (zigzag 
ribbon) exhibits unique electronic states with energies 
close to the Fermi level. These states are highly local- 
ized at the edge carbon atoms, and are defined as edge 
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states, but the states do not appear in the armchair 
ribbonsi 18 ' 19 i 20 i 21 i 22 Electronic^ magnetic; 18 i 19 i 24 opti- 
cal^ and field emissionSS»21 properties have also been 
studied extensively with respect to potential edge ef- 
fects. However, to the best of the authors' knowledge, 
the phonon transport properties of graphitic ribbon at 
low temperature have not been studied. Despite a long- 
standing theoretical interest going back to Peierls' early 
workj 2 ^ little progress has been made toward elucidat- 
ing phonon transport in mesoscopic or nanoscale systems 
with dimensions much smaller than the mean free path 
of phonons. Recent advances in nanotcchnology have 
made it possible to measure the thermal conductance 
by phonons in such small systems and to observe the 
universal quantum of thermal conductance. The value 
obtained, n 2 k 2 B T/3h = 9.4 x 1(T 13 [W/K]^S is consis- 
tent with the value predicted theoretically by Rego and 
Kirczenow^i and by other theoretical studies i22i2ii Very 
recently, we reported that the quantization of thermal 
conductance in SWNTs is independent of tube radius or 
chirality^ Motivated by these discoveries, the objective 
of the present study is to calculate the phonon dispersion 
relations based on the empirical force-constant model and 
to demonstrate the universal features of quantized ther- 
mal conductance in graphitic ribbon at low temperature. 

This paper is organized as follows. In section [HI 
the force-constant method based on Brenner's empirical 
potential is introduced for investigation of the thermal 
properties of carbon nanostructures, focusing on SWNTs 
and graphitic ribbons as primary examples. To ver- 
ify the transferability of the method, the phonon dis- 
persion relations and the density of states of SWNTs 
are calculated. In section IIHI the atomic structures of 
graphitic ribbon are optimized using the Brenner poten- 
tial and the phonon dispersion relations of graphitic rib- 
bon are obtained by the empirical force-constant method. 
In section IIVI phonon-derived thermal conductance in 
graphitic ribbon is discussed based on theoretical analy- 
sis using the Landauer theory of phonon transport, and 
the quantization and universal behavior of thermal con- 
ductance in graphitic ribbon are presented. The paper is 
finally summarized in section Ivl 



II. EMPIRICAL FORCE-CONSTANT MODEL 

The phonon dispersion relations of a system are ob- 
tained by diagonalizing the dynamical matrix given by 

D£{k) = £ — l==^exp (ife • (J^ - R n )) . (1) 

Here, the summation index m runs over all unit cells, k is 
the wave vector, R m is the spatial coordinate of the mth 
unit cell, and Mj is the mass of the ith atom. D^(k) is 
independent of n owing to the translational symmetry of 
the system. The force-constant is defined as the 

second derivative of the potential with respect to atomic 



coordinates as follows. 

FP-V 

K mjl3 = — (2) 

ma Br a Br -a' 

where V is the potential energy of the system, r n i a is the 
a component (a = x, y, z) of the position vector of the 
ith atom in a unit cell n. 

For carbon nanostructures, Brenner's empirical poten- 
tial is adopted for V in Eq. (|2J). The Brenner potential 
has been successfully used for MD simulations of carbon 
nanostructures. Details of this potential can be found 
in Brenner's original paper^ The force constants can of 
course be obtained for small systems from the Hellmann- 
Feynmann force computed by first-principles calculations 
based on density functional theory. However, for large 
systems, the phonon transport properties must be cal- 
culated for systems with a very large phonon mean-free- 
path, for example, the phonon mean-free-path of SWNTs 
is of the order of 0.1 — 1.0 /im^ The empirical method, 
on the other hand, can be computed very quickly, mak- 
ing it very useful for computing atomistic quantities in 
larger systems. 

Substituting the Brenner potential into Eq. J2J , an an- 
alytical form of the force-constant can be obtained in 
a straightforward manner. The derivation is omitted 
here for the sake of brevity. The reliability and useful- 
ness of the analytical equation for analysis of the ther- 
mal properties of carbon nanostructures is evaluated by 
comparing phonon dispersion relations for SWNTs de- 
termined by this model with the results present in previ- 
ous studies. As the force constants are only determined 
for the equilibrium positions of carbon atoms in the sys- 
tem, it is necessary to optimize the atomic geometry of 
the SWNTs using the Brenner potential. For an SWNT 
with a chiral vector of Ch — (10,10), where the chiral 
vector Ch = (n, m) uniquely determines the geometri- 
cal structure of the SWNT ) 36 i 37 the optimal bond length 
between two carbon atoms is obtained as 1.45 A using 
the parameter set given in Table III in Brenner's original 
paper 3 . This is slightly longer than the typical exper- 
imental value of 1.44 A. Figure ^a) shows the phonon 
dispersion relations for the (10, 10) SWNT, which has 40 
carbon atoms in the unit cell and 120 vibrational degrees 
of freedom. 

Due to mode degeneracies, only 66 distinct phonon 
branches including 12 nondegenerate branches and 54 
doubly degenerate branches appear in Fig. da). This 
result is consistent with the number of distinct phonon 
branches obtained by point group theory for atoms in the 
unit cell. The corresponding density of states is shown 
in Fig.^b). The small peaks represent van Hove singu- 
larities unique to one-dimensional systems. 

Four acoustic modes with linear dispersion appear in 
the low-frequency region of Fig. ^a): a longitudinal 
acoustic (LA) mode, a twisting acoustic (TW) mode, 
doubly degenerate transverse acoustic (TA) modes. The 
sound velocities of the LA, TW, and TA modes for 
(10, 10)-SWNT are estimated to be 22.2, 13.3, and 
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10.3 km/s, respectively. The energy gap of the doubly 
degenerated lowest optical modes depends only on the 
tube radius R, and decreases approximately according to 
~ 1/R 2 . These results are in good agreement with pre- 
viously reported phonon dispersion relations for SWNTs 
determined by a dynamical matrix method with force- 
constant parameters for graphene^ and first-principles 
calculations Therefore, the phonon dispersion rela- 
tions for other carbon nanostrucures are considered to be 
reliably obtained using the force-constant method based 
on Brenner's empirical potential. 



III. PHONON DISPERSION RELATIONS FOR 
GRAPHITIC RIBBON 

The empirical force-constant method is employed to 
calculate the phonon dispersion relations of the two types 
of graphitic ribbon without hydrogen termination. The 
dehydrogenated armchair- and zigzag-ribbons are sta- 
ble even at high temperatures above 2000K. 38 Figure [3 
shows the schematic structures of the dehydrogenated 
armchair- and zigzag- ribbons. Here, N a and N z denote 
the number of dimcr lines and zigzag lines, respectively. 
Before calculating the phonon dispersion relations, the 
atomic structures need to be fully optimized using the 
Brenner potential. For the armchair ribbon with N a = 9 
(18 atoms per unit cell), the optimized bond length be- 
tween two edge atoms is 1.39 A, which is shorter than the 
typical bond length of 1.45 A inside the armchair ribbon. 
This shortening of the edge bond length appeared for 
all armchair ribbons calculated in this study, and results 
from the fact that two atoms at the edges of the arm- 
chair ribbon form a triple bond, C=C. This triple-bond 
feature has been already revealed in the first-principles 
calculations of the electron density of armchair ribbons 
without hydrogen termination^ 

The phonon dispersion relations represent key informa- 
tion for identifying the triple bond at the edge of the arm- 
chair ribbons experimentally. As a result of the strong 
triple bond of C=C, doubly degenerated branches with 
frequency of approximately 1800 cm -1 at the T point 
appear in the phonon dispersion relations. These two in- 
plane modes correspond to antiphase and in-phase C=C 
stretching motions on both edges of the ribbon. The fre- 
quency of ~1800 cm" 1 is higher than the typical SWNT 
peak at ^1590 cm" 1 in the Raman spectrum originating 
from the Raman- active Ei g mode of a graphitic sheets 
As this is indicative of a Raman shift at this frequency, 
it is possible to identify the C=C bond at the edge of an 
armchair ribbon based on the Raman spectrum. 

The edge-localized phonon modes appearing as out-of- 
plane modes of the armchair ribbon are shown in Fig. 
for three different ribbon widths, together with the dis- 
persion relations for a graphene sheet projected onto the 
armchair axis (parallel to the translation vector along the 
ribbon axis a) . The upper and lower shaded areas corre- 
spond to the optical and acoustic modes of the graphene 



sheet, respectively. For all the armchair ribbons calcu- 
lated, four special branches were found outside the region 
of the projected dispersion relations of a graphene sheet, 
as indicated by arrows in the figure. The existence of 
these branches has been already predicted by Igami et 
ul.^iLlL but the detailed properties have yet to be re- 
ported. From an examination of the atomic displace- 
ments, the phonon states in these four branches appear 
to correspond to a strongly localized state on the arm- 
chair edges. Figure 01 shows snapshots of the atomic dis- 
placements at k\a\/ir = 1 for these edge-localized phonon 
states of the armchair ribbon with N a = 9. As shown 
in Fig. |3| two branches occur below the lower shaded 
area, representing an acoustic modes and a lowest optical 
mode. In the acoustic mode, both edge-dimers in the unit 
cell undergo in-phase vibration (Fig. Ufa)), while in the 
lowest optical mode, both edge-dimers exhibit anti-phase 
vibration (Fig. |I£b)). Two other edge-localized modes 
appear in the region between the upper and lower shaded 
areas in Fig. |3| representing a mode due to two carbon 
atoms in the anti-phase vibrating edge-dimer and the in- 
phase vibration of the two edge-dimers on each side in 
the unit cell, as shown in Fig.^Jc), and an edge-localized 
mode in which both edge-dimers vibrate in anti-phase, 
as shown in Fig. 0Jd). 

Figure [5] shows the phonon dispersion relations of the 
out-of-plane modes for three zigzag ribbons, together 
with the dispersion relations (shaded areas) of a graphene 
sheet projected onto the zigzag axis (parallel to the vec- 
tor z). The upper and lower shaded areas contact at 
k\z\/n = 2/3. The two branches appearing immediately 
below the lower shaded area cannot be predicted from 
the projected dispersion relations of a graphene sheet. 
The phonon states in these two branches correspond to 
edge-localized states from observations of the atomic mo- 
tions in these branches (Fig. EJ). In the acoustic mode, 
two edge-atoms in the unit cell exhibit in-phase vibra- 
tion (Fig. |Ua)), while in the lowest optical mode, the 
two edge-atoms vibrate in anti-phase (Fig.HJb)). 

The lowest-lying phonon dispersion relations govern- 
ing the low-temperature thermal-transport properties of 
graphitic ribbon are represented by three acoustic modes 
in the long-wavelength limit; an LA mode, an in-plane 
TA mode, and an out-of-plane TA mode. In contrast to 
SWNTs, the two TA modes are not degenerate because 
there is no rotational symmetry around the LA-wave di- 
rection. Using the empirical force-constant method, the 
energy gap of the lowest optical mode at the T point, 
A_E op = fia> p, was shown to depend only on the rib- 
bon width W according to the relation AE op oc 1/W. 
That is, the energy gap is independent of the edge shape. 
The AE op of SWNTs has an approximately quadric ra- 
dius dependence given by 1/R 2 for tube radius in 
contrast to the AE op oc 1/R relation predicted by the 
zone-folding approach^ This difference arises because 
the zone-folding approach neglects the effects of cylinder 
geometry or curvature of the SWNT. AE op represents 
an important physical parameter for characterizing the 
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quantization plateau of thermal conductance, as shown 
in the next section. 



IV. PHONON TRANSPORT IN GRAPHITIC 
RIBBON 

From an experimental fact that the phonon mean-free- 
path of SWNT is 0.5-1.5/xm in the wide temperature 
range from 350K to 8K^ we can expect the graphitic rib- 
bon is a ballistic phonon conductor with a large phonon 
mean-free-path of the order of 1/im at low temperatures. 
Consider a situation in which phonons travel ballistically 
through a graphitic ribbon, connected smoothly between 
two thermal reservoirs with a small temperature differ- 
ence. 

It has been shown previously in a study of the thermal 
conductance of SWNTs 34 that the thermal conductance 
in such situation is well described by 

9k 2 T , 

<T) = -f- £ n ) - , (3) 

m 

based on the Landauer theory of phonon transport. 
Here, G(x) is defined as 

G(x) = 0(2, e-«) + x4>(l, e- x ) + y /(*), (4) 

i min/max . i . i r 

where x m is the minimum/ maximum value ot x m = 
hu m /kBT, <fi(z, s) — X)^Li( s ™/ rt2 ) i s the Appel function, 
and f(x) — l/(e x — 1) is the Bose-Einstein distribution 
function. Details of the derivation of Eq. © can be found 
in the previous report^ 

Substituting w™ m and w™ ax for the graphitic ribbon 
into Eq. @ gives the thermal conductance of these two 
types of graphitic ribbon. FigureEfa) shows the temper- 
ature dependence of the thermal conductance of graphitic 
ribbon with different widths and edge shapes, normalized 
against a universal value of 3kq (k = ir 2 k 2 B T /3h is the 
universal quantum of thermal conductance). All calcu- 
lated curves approach unity in the low-temperature limit, 
indicating that the thermal conductance in graphitic rib- 
bon is quantized at the universal value of 3kq, indepen- 
dent of the ribbon width or edge shape. It can be seen 
from Eq. © that each acoustic mode with x™ m = 
(LA mode and in-plane and out-of-plane TA modes) con- 
tributes a universal quantum of kq to the thermal con- 
ductance. 

The quantization plateau seen in Fig. Ela) is broken 
with increasing temperature, with the thermal conduc- 
tance beginning to rise as the optical phonons become 
excited. Since the characteristic energy scale of the 
graphitic ribbon in this temperature region is the en- 
ergy gap of the lowest optical mode, AE op , it would 
be convenient to introduce a dimensionless temperature 
t = fcgT/ AE op . Taking account of the three acoustic 
modes and the lowest optical mode and substituting the 



value of w™ m for these branches into Eq. @ , the thermal 
conductance can be given in the following simple form. 

This equation gives the thermal conductance as being 
dependent only on r. As shown in Fig. 17Tb). all curves 
in Fig. 0a) collapse onto a single curve of Eq. JSJ once 
plotted against r. The single curve rises at r rj 0.1 
from the quantization plateau. The width of the quan- 
tization plateau with respect to temperature, AT w 
0.1 x AE p/kB, is found to be inversely proportional to 
the ribbon width W, that is, AT cx 1/W. This means 
that the phonon transport properties of graphitic ribbon 
at low temperature are independent of the atomic struc- 
ture such as edge shape. This follows from that fact that 
the wavelength of phonons as heat carriers is sufficiently 
high in comparison with the bond length of the graphitic 
ribbon. 

Inspired by recent experiments on phonon transport in 
SWNTsjS we previously showed that the phonon-derived 
thermal conductance of SWNTs at low temperatures is 
quantized at 4kq, corresponding to 4 acoustic modes; the 
LA mode, two TA modes, and the TW modei^i Since 
there is no TW mode in the phonon dispersion curves of 
graphitic ribbon, the thermal conductance in graphitic 
ribbon is smaller than that in SWNTs. In graphitic rib- 
bon, one of the three quanta is carried by edge-localized 
acoustic phonons. Furthermore, while the quantization 
plateau width for graphitic ribbon extends in proportion 
to 1/W with decreasing W, independent of edge struc- 
tures, that for SWNTs is proportional to 1/R 2 . There- 
fore, the width of the quantization plateau for narrow 
graphitic ribbon is smaller than that for SWNTs with 
the same number of atoms per unit cell. 

V. SUMMARY 

As a first step toward systematic study of the ther- 
mal properties of carbon nanostructures, the phonon 
dispersion relations and density of states of carbon 
nanostructures were calculated using a force-constant 
method based on Brenner's empirical potential, focus- 
ing on SWNTs as a primary example. The results ob- 
tained are in good agreement with previous results for 
the phonon dispersion relations of SWNTs determined 
by first-principles calculations ^iLii The empirical force- 
constant method is therefore considered to be a reli- 
able method for calculations of the phonon dispersion 
relations for various carbon nanostructures. Using this 
method, phonon dispersion relations were calculated for 
graphitic ribbon with armchair- or zigzag-shaped edges, 
and localized phonon modes were identified at the edges 
of the ribbons. For armchair ribbons, two in-planc 
phonon modes with high frequency (near 1800 cm -1 ) 
are assigned to C=C stretching motion at the edges. 
The out-of-plane TA and lowest optical modes were also 
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shown to be edge-localized phonon states in the arm- 
chair and zigzag ribbons. These low- lying modes were 
found to contribute to the thermal transport properties 
of graphitic ribbon at low temperature. 

Using the Landauer formula for phonon transport, per- 
fect transmission of an acoustic mode in graphitic rib- 
bon was found to contribute a universal quantum of 
kq = Tr 2 kgT/3h to thermal conductance. The width 
of the temperature-dependent quantization plateau was 
also observed to be inversely proportional to ribbon 
width. Quantization phenomena of thermal conductance 
are therefore expected to be experimentally observable 
in various one-dimensional carbon nanostructures such 
as carbynes and polyacetylenes as well as SWNTs and 
graphitic ribbons. The successful application of an em- 



pirical force-constant method combined with the Lan- 
dauer formula in this study revealed the quantized and 
universal features of phonon-derived thermal conduc- 
tance in graphitic ribbon, demonstrating the essential 
role of ballistic edge-phonon modes in these materials. 
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FIG. 1: (a) Phonon dispersion curves and (b) density of 
states for a single-walled carbon nanotube with chiral vector 
(10, 10). \T\ denotes the magnitude of the translation vector 
along the tube axis. 

FIG. 2: Atomic structure of (a) armchair ribbon (N a = 7) 
and (b) zigzag ribbon (N z = 4). Broken rectangle represents 
a unit cell, and a and z denote the translation vectors along 
the ribbon axis. 



FIG. 3: Phonon dispersion relations of armchair ribbon with 
N a = 9, 11, and 15. Shaded areas represent the phonon 
dispersion relations of a two-dimensional graphitic sheet pro- 
jected onto the armchair axis. Edge-localized phonon modes 
are indicated by arrows. 



FIG. 4: Atomic displacement of edge-localized phonon states 
at k\a\/ir — 1 in the out-of-plane modes of an armchair rib- 
bon of N a = 9 for (a) the transverse acoustic mode, and (b) 
the lowest-lying optical mode, (c), (d) Atomic displacement 
of states in the two branches between the upper and lower 
shaded areas in Fig. [3] Solid circles represent edge carbon 
atoms. 



FIG. 5: Phonon dispersion relations of zigzag ribbons with 
N a = 8, 10, and 16. Shaded areas represent the phonon 
dispersion relations of a two-dimensional graphitic sheet pro- 
jected onto the zigzag axis. Edge-localized phonon modes are 
indicated by arrows. 



FIG. 6: Atomic displacements of edge-localized phonon states 
at k\z\/n = 1 in the out-of-plane modes of zigzag ribbon with 
N a = 6 for (a) the transverse acoustic mode and (b) the 
lowest-lying optical mode. Solid circles represent edge carbon 
atoms. 



FIG. 7: (a) Low-temperature phonon-derived thermal con- 
ductance in graphitic ribbon. Open circles denote armchair 
ribbon (N a ~ 9, 11, and 15), and solid circles denote zigzag 
ribbon (N z — 8, 10, and 16). (b) Thermal conductances of 
armchair (N a = 11 and 15) and zigzag (N a = 10 and 16) rib- 
bon as a function of temperature scaled by the energy gap of 
the lowest optical mode. Broken curve represents a universal 
curve given by Eq. @. 
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